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ABSTRACT 

We discuss the present state of knowledge and thought concerning the spread in age found 
among Galactic globular clusters, with some discussion of the implications for what hap- 
pened during the earliest stages of the formation of the Milky Way Galaxy. Differential 
observational techniques to derive the relative age differences among clusters of similar 
metallicity are discussed in detail. We conclude that as of the current date (mid-April 
1996) the state of the field is still somewhat muddled. However, we believe that there is 
now a substantial body of evidence — including a particularly revealing intercomparison 
of the color-magnitude diagrams of NGC 1851, NGC288, and NGC362 presented here - 
indicating that age is not the dominant second parameter determining the shape of glob- 
ular clusters' horizontal branches. If our assertion is correct, then apart from a handful of 
anomalous clusters that may well have been captured from a satellite dwarf galaxy, there 
is no strong evidence either for a significant spread in age among clusters of a given metal 
abundance or for a systematic variation of mean age with Galactocentric distance. On the 
question of whether there is a significant age difference between metal-poor and metal-rich 
clusters, we feel compelled to fall back on the Scottish verdict: "Unproven." Data now 
being collected by numerous groups in various subdisciplines may resolve the remaining 
controversy within a few years. 

1. INTRODUCTION 

Few numbers in science are more avidly sought than the ages of the Galactic globular 
clusters (GCs). These spherical, compact, gravitationally bound systems of typically 10 5 - 
10 6 stars that orbit the center of the Milky Way (and most, if not all, other large galaxies), 
are believed to be among the first recognizeable stellar systems that formed in the Universe 
(see Peebles and Dicke 1968). As they are the most ancient objects known for which reliable 
ages can be determined, and as the Universe cannot be younger than the oldest object it 
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contains, the GCs provide one of the most important of the few constraints that we have 
on cosmological models. In tandem with this study, VandenBerg, Bolte, and Stetson 
(1996) review our present understanding of absolute cluster ages and briefly assess their 
implications for standard Friedmann cosmologies. 

The calibration of absolute ages is still subject to observational and theoretical un- 
certainties at the ~ 20% level, and represents a major limitation on our ability to test 
cosmological models. Nevertheless, it is quite possible to determine relative GC ages with 
sufficient precision to address some outstanding questions concerning the formation of our 
own Milky Way Galaxy. In particular, we now have the means to estimate — to within a 
fraction of a Gyr — the dispersion in age of globular clusters having very similar chemical 
compositions, as well as the variation of age with Galactocentric distance encompassed 
by (once again) those systems of comparable metal abundance. Even though the deriva- 
tion of a reliable age-metallicity relation for the GCs remains problematic, due largely to 
distance-scale uncertainties, we may soon be able to ascertain whether the protoGalaxy 
underwent the type of rapid, monolithic collapse envisioned by Eggen, Lynden-Bell, and 
Sandage (1962), or whether the Galaxy formed much more gradually through the accretion 
of independent sub-systems over a prolonged period of time (cf. Searle and Zinn 1978). 
It is the purpose of this investigation to review the many advances that have been made 
in recent years towards measuring relative ages of GCs with a precision sufficiently high 
(^ lGyr) to resolve the first epoch of star formation in the Galaxy. Of course, there 
are many other paths of investigation that are leading toward a complete understanding 
of our Galaxy's childhood and adolescence based on studies of halo and thick disk field 
stars; investigations of the Sagittarius galaxy, Carina, and other nearby dwarfs; and in situ 
observations of galaxies "now" forming at large redshift. While we will not discuss these 
approaches in this review, we do not in any way wish to minimize their importance. 
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2. BACKGROUND: CONCEPTS OF HOW THE GALAXY FORMED 
The general notion that the Galaxy reached its present form after a self-gravity-driven 
collapse from a more dispersed state dates back at least to Kant (1755). This classical view 
of the formation of the Galaxy has been given a modern codification and quantification by 
Eggen, Lynden-Bell, and Sandage (1962, hereafter ELS), who proposed that the original 
protoGalaxy consisted of a single large gas cloud with a radius (in present-day coordinates) 
of order 100 kpc. Unable to support itself by gas pressure once its mass had stopped the 
local Hubble expansion, this cloud immediately underwent a free-fall collapse. The globular 
clusters and field halo stars are presumed to have formed in condensations embedded within 
the general intraGalactic medium during this rapid free-fall phase, whose duration could 
not have exceeded a few times 10 8 years. Finally, when the gas achieved a density which 
allowed pressure support, and dissipation therefore began to become important, it was 
able to radiate away its kinetic energy; residual angular momentum then forced it to settle 
into a rotating disk, where all subsequent star and cluster formation has taken place. In 
contrast to the gas, the dissipationless halo field stars and GCs have retained the kinetic 
energies, angular momenta, and general spatial distribution that they possessed at the 
instant the gaseous collapse ceased. 

The principal alternative to this model was presented by Searle and Zinn (1978, here- 
after SZ). Arguing from the postulate that differences in age were responsible for the range 
in horizontal-branch types at fixed metal abundance — the so-called "second parameter" 
effect (van den Bergh 1965, 1967; Faulkner 1966; Sandage and Wildey 1967) — they in- 
ferred that the formation of the Galactic halo could not have been as rapid as a free-fall 
collapse. According to early horizontal- branch models (e.g., Faulkner 1966; Rood and 
Iben 1968; Simoda and Iben 1970; Castellani and Tornambe 1977), the observed range 
in horizontal-branch morphological types would require that the epoch of globular-cluster 
formation must have been 3>1 Gyr. SZ achieved the required prolonging of the cluster 
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formation time in the outer halo by abandoning the notion of a continuous medium in 
free-fall radial collapse. Instead, they adduced the theoretical models of Binney (1976), 
which assumed that the protoGalactic material consisted of distinct "flows" of gaseous 
material which would not be directed purely in the radial direction, but which would on 
occasion collide with each other, dissipating their kinetic energy by shock-induced radia- 
tion and producing transient sheets of cool, compressed gas suitable for the formation of 
a generation of stars and clusters. If the filling factor of these gas flows was sufficiently 
small, the time between collisions could be long enough to spread out the dissipation of 
the original kinetic energy over the time span required to produce the inferred age spread 
in the outer halo. 

The relationship between these transient sheets forming as a result of collisions be- 
tween currents of protoGalactic gas, and the "protogalactic fragments" discussed earlier 
by Searle (1975), is not entirely clear. Searle's earlier paper expands upon previous discus- 
sions by von Weizsacker (1955) and Oort (1958), who may have been the first to suggest 
that the collapse of the Galaxy was not the simple radial, gravity-driven collapse envi- 
sioned by Kant and by ELS, but rather that turbulent motions within the protoGalactic 
material moderated the rate of contraction through the gradual dissipation of kinetic en- 
ergy. Searle originally adopted the speculations of von Weizsacker and Oort to explain 
why the range of metal abundance among stars within a given GC is immeasureably small 
compared with the total chemical enrichment in the typical cluster. This circumstance 
would be extremely hard to understand if the material which became a given globular 
cluster separated itself out from the general protoGalactic medium early in the free-fall 
phase, and then enriched itself in heavy elements through successive supernova events. One 
would then expect to see a range of metal abundances from zero up to some maximum, 
representing stars that formed in pockets of variously enriched gas during the intervals 
between the supernova explosions. Instead, Searle suggested that the protohalo consisted 
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of "initially gaseous fragments, each small compared with the Galaxy but large compared 
with globular clusters. . . . Following Kaufman's (1975) discussion of the efficacy of mixing 
by supernova explosions, it seems reasonable to assume complete mixing within a fragment 
and no mixing between them." Thus, within a given fragment, initial generations of stars 
would produce chemical enrichment and complete remixing of those secondary elements 
throughout the mass of the fragment. At various times, star clusters would form, each 
inheriting the unique metallicity belonging to the parent fragment at that instant. Each 
fragment would eventually encounter the main body of the Galaxy; its gaseous content 
would be stripped and added to the Galactic disk, while its stars and clusters would return 
to the halo. Subsequently, Searle and Zinn proposed the very similar (though not identi- 
cal) scenario already discussed to enable them to adopt a very different postulate: namely, 
that the second parameter is age. The similarity of the two situations envisioned by Searle 
and by SZ has led some authors — including, on occasion, us — to confuse the two, to 
speak of "Searle-Zinn fragments" and to suggest that the present-day dwarf spheroidal 
galaxies may be the aged remains of a few of these primordial lumps of protoGalactic 
material that have managed to survive to the present day. So does the model really en- 
vision many independent long-lived condensations zipping through otherwise empty space 
like bees buzzing around a hive, all the while undergoing continuous internal evolution 
until, one by one, they collide and merge their gaseous bodies to form the kinematically 
ordered Galactic disk, while their collisonless stellar and cluster progeny continue to roam 
the original spherical volume once occupied by the whole protoGalaxy? Or were the stars 
and clusters formed in transient sheets of cool material behind the shock fronts of colliding 
gas streams, sheets which formed but one generation of stars and clusters before they were 
ionized by OB stars and torn apart by supernovae and returned to the turbulent, gaseous, 
spheroidal protoGalaxy — a process that has nothing whatever to do with modern dwarf 
galaxies? Even Oort spoke alternatively of "irregular streamings" and of "gas clouds" 
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which undergo collisions; it is not clear that Searle and Zinn specify which version they 
prefer, either. 

More recently, the detailed theoretical modelling of horizontal branch morphologies 
by Lee, Demarque, and Zinn (1988, 1994; also Lee 1991, 1992, 1993) have suggested that 
the time span required for the formation of halo objects is ^,5 Gyr; in particular, Lee, 
Demarque, and Zinn (1988) inferred an age difference of 6 Gyr for the well-known second- 
parameter pair of globular clusters, NGC 288 and NGC362, if the absolute age of the 
older of the two (NGC 288) is 18.5 Gyr — a difference-over-average of ~40%. In their 
more recent discussion, Lee, Demarque, and Zinn (1994) revised the estimated difference 
to 3-4 Gyr for an adopted age of 14.9 Gyr for NGC 288 — a difference of ~30%. 

Sandage (1990) has argued for what is essentially a synthesis of the two models. He 
contends that critics of ELS have been overly literal in their interpretation of its highly 
simplified picture, as naturally there would necessarily have been some variation of local 
densities within the protoGalactic gas cloud, and denser lumps would have collapsed in 
upon themselves on shorter timescales even as the Galaxy overall collapsed in upon itself on 
a longer timescale. Still, there remains a dichotomy of opinion, with some agreeing with 
Sandage that the collapse was "coherent" (albeit noisy), "rapid" (in the sense that the 
radial infall greatly exceeded turbulent motions, as required by the highly eccentric orbits 
of the halo stars), and dominated by gravity; while others agree with SZ that, at least in 
the outer halo, the collapse was "chaotic" and "slow" (required by the age interpretation 
of the second-parameter problem), and that gravity initially dominated over kinetic energy 
by only a tiny margin. The most extreme version of the latter picture is probably that 
of Zinn (1993), who argued that the "young" globular clusters in the outer halo had their 
origin in just a few satellite galaxies (resembling, perhaps, the Magellanic Clouds) which 
have only quite recently been disrupted and accreted by the Galaxy. 
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It must be recalled that, up to 1987, the age explanation of the second-parameter 
phenomenon was still largely a hypothesis, relying upon assumptions concerning the uni- 
versality of chemical abundance ratios, the internal angular momentum of stars, mass-loss 
mechanisms, and the like. None of the so-called "second-parameter" clusters had had its 
chronological age accurately measured by the more fundamental techniques based on the 
luminosity and temperature of the main-sequence turnoff. But, beginning in 1987, such 
measurements began to appear in rapid order. 

3. MODERN MEASUREMENTS OF RELATIVE CLUSTER AGES 
The true measure of a star cluster's age is the intrinsic luminosity of its main-sequence 
turnoff. A star can be thought of as a gravitationally confined reservoir of nuclear fuel. 
During nearly all of its lifetime, the luminosity of a star is set by the requirement that it 
support itself against gravitational contraction, and a new, young star quickly adjusts 
its internal density and temperature structure until the fusion of hydrogen to helium 
in the stellar core provides precisely the required rate of energy production. Over the 
approximately 90% of a star's total lifetime that is spent in this so-called "main-sequence" 
phase of its evolution, it grows gradually more luminous and slightly hotter in response 
to the changing chemical structure of its interior. There comes a time in the life of a 
star, however, when it has exhausted the available hydrogen fuel at its center. Then 
the inert helium core, containing some 10% of the star's mass, shrinks to a volume of 
order the dimensions of the planet Earth while the outer layers of the star rapidly expand 
in radius, causing the surface temperature of the star to fall precipitously as hydrogen 
continues to be consumed in an energy-generating shell around the helium core. The 
rate of fuel consumption during the main-sequence phase of evolution is known from the 
star's luminosity; the total amount of fuel available in the hydrogen-burning core is easily 
estimated from even the most primitive stellar interior models; and the ratio of these two 
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numbers yields the amount of time a star of a given mass can remain on the main sequence. 
More massive stars consume their fuel far more rapidly than less massive stars, and die 
proportionately sooner. A star cluster represents a coeval assemblage of stars spanning a 
range of mass: by simply observing a cluster and determining the mass or luminosity of 
stars that are just switching over from slow hotward evolution to rapid coolward evolution, 
one has a direct measure of the age of the cluster. 

Although the commonly held wisdom through the mid-1970's was that globular clus- 
ters were coeval, this interpretation was based on photoelectrically calibrated photographic 
photometry. The difficulties of photometering main-sequence-turnoff and fainter stars on 
photographic plastes were formidable. In retrospect it is clear that even age dfferences as 
large as 3 - 4 Gyr could have been lost in the random and systematic photometry errors. 
However, when CCDs came into common use in the mid-1980s and the software for de- 
convolving overlapping stellar images in the invariably crowded GC fields became widely 
available, the stage was set for establishing relative GC ages with uncertainties < 1 Gyr. 

The first clear measurement of an anomalously young age for a Galactic globular clus- 
ter emerged from two studies of the outer-halo, second-parameter cluster Palomar 12, which 
Ortolani (1987) reported to be younger by some substantial (but unspecified) amount than 
the other Galactic GCs; at the same conference, Stetson and Smith (1987) independently 
reported an anomalously young age for Pal 12, and estimated by a differential comparison 
with 47 Tucanae and M 5 that it was some 30% younger than those two nominally normal 
clusters. Subsequent more exact analyses (Gratton and Ortolani 1988; Stetson et al. 1988, 
1989) supported this provisional estimate. Since then, the clusters Ruprecht 106, Terzan 7, 
Arp 2, and IC 4499 (Buonanno et al. 1990, 1993, 1994, 1995; Kubiak 1991; Da Costa, 
Armandroff, and Norris 1992; Ferraro et al. 1995) have been added to the roster of clusters 
unequivocally younger than the majority. However, as discussed by Lin and Richer (1992) 
and Fusi Pecci et al. (1995), it is possible that these clusters have been stripped from the 
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Magellanic Clouds by the tidal force of the Galaxy, or from another satellite galaxy that 
has since been disrupted. Indeed, the recent discovery of the dwarf spheroidal galaxy in 
Sagittarius (Ibata, Gilmore, and Irwin 1994), and the possibility that the globular clusters 
NGC 6715 = M54, Arp 2, Terzan 8, and conceivably Terzan 7, may be associated with it, 
could be a contemporary example of this process. Therefore it is unclear whether these 
clusters trace the dominant formation mode of the original halo, or just a small minority 
component of it. 

To understand the principal formation mode of the Galactic halo, it will be necessary 
to measure precise relative ages for a majority of the globular clusters. Work to date tends 
to suggest that the age dispersion among the nearest clusters is small, but perhaps per- 
ceptible. Bolte (1989) obtained precise relative photometry for the main-sequence turnoffs 
of the well-known nearby second-parameter pair NGC 288 and NGC 362, from which he 
concluded that the former was older than the latter by some 3 Gyr. Green and Norris 

(1990) confirmed this result shortly thereafter (however, see §4 below). Note that this age 
difference is appreciably smaller than the one originally suggested by Lee et al. (1988): 
they had inferred an age of 18.5 Gyr for NGC 288, versus 11.2 Gyr for NGC 362. Lee 

(1991) subsequently pointed out that his previous estimate of the age difference had been 
made on the basis of a Solar ratio of oxygen to iron in the globular-cluster stars, and that 
allowing for various effects related to an enhancement of oxygen would reduce the absolute 
size of the age difference. Their most recent calibration of theoretical horizontal-branch 
morphologies (Lee et al. 1994) yields a relative age difference of 3-4 Gyr, if the absolute 
age of NGC 288 is 14.9 Gyr. This would appear to be a firm lower limit to the difference 
in age between these two clusters if age is the sole cause of the disparity in their respective 
HB morphologies. Catelan and de Freitas Pacheco (1993, 1994) have found from their HB 
simulations that an age difference as small as w 2-3 Gyr would require that either both 
clusters are ^ 10 Gyr old or some parameter in addition to age is different between them. 
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VandenBerg, Bolte, and Stetson (1990; hereafter VBS) compared the color-magnitude 
diagrams ("CMDs") of another nearby second-parameter pair, NGC 5272 = M3 and 
NGC 6205 = M 13, and estimated an age difference of some 1-2 Gyr (7-15%) as com- 
pared to the Lee et al. (1994) prediction of 3 Gyr (>20%) from their synthetic HB 
calculations. Once again, it would seem to be impossible for modern HB models to accom- 
modate a smaller difference in age between M 3 and M 13 — should that turn out to be 
the case — unless something other than, or in addition to, age is varying (see Catelan and 
de Freitas Pacheco 1995). On the other hand, VBS found no detectable age range in excess 
of some 0.3 Gyr (~2%) among the most metal-poor clusters with the best photometry; 
these particular results thus tend to favor ELS, especially the Sandage (1990) version. 

Therefore, the evidence from nearby clusters suggests that understanding the forma- 
tion of the Galactic halo is not so much a matter of deciding which of ELS and SZ is 
"right," but rather accepting that there is some truth in both pictures, and trying to de- 
cide the fraction of the halo that is best described by each of them. About six orders of 
magnitude separate the mass of the typical globular cluster from the mass of the Galaxy. 
Where within this range did the mass of the typical cluster-forming lump lie? What was 
the ratio of kinetic to gravitational potential energy in the protoGalaxy at the instant the 
local Hubble expansion ceased and the contraction began? Before we undertake to answer 
these questions, it must be stressed that the observational studies discussed so far have 
had little to say about clusters far from the Solar circle, and thus are almost certainly 
unfairly sampling the full range of circumstances that prevailed in the protoGalactic halo. 
To obtain the big picture, it will be necessary to measure precise relative ages for clusters 
sampling different extremes of ancestry: those close in to the Galactic center, and those at 
the outermost fringes of the modern-day halo. Reddening, crowding, and confusion with 
field stars complicate observations of the former, while sheer distance makes the latter 
difficult to study in any detail. But, before progressing to a discussion of the most recent 
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age determinations for globular clusters, it is worthwhile making a slight diversion to a 
discussion of how relative ages can be measured. 

4. DIFFERENTIAL AGE-DATING TECHNIQUES 
Since the two dominant competing models for the formation of the Galaxy make 
conflicting predictions for the relative duration of the halo-formation epoch (i.e., ~1% 
in classical ELS or ^10% in Sandage's (1990) more recent reformulation, as compared 
to ~30% in classical SZ and 30-100% in Zinn's (1993) most recent extension of that 
model), precise relative age determinations should suffice to distinguish between the two 
extremes of possibility. This is fortunate, because relative age determinations can use 
stellar evolution theory in a strictly differential sense, removing most of the effects of 
theoretical uncertainties in absolute chemical abundance ratios, opacities, convection for- 
malism, temperature-color relations, and the like. Differential comparisons can also be 
devised which reduce the effects of observational uncertainties in the absolute distance 
scale, overall metal abundance, and individual cluster reddenings. 

The differential age-dating techniques for globular clusters in use at the present time 
can be separated into two basic classes: those that are "vertical" and those that are 
"horizontal" in the color-magnitude diagram. Each has its advantages and its problems. 
The most venerable of the vertical techniques is the magnitude difference between the 
luminosity of the horizontal branch — as measured from the RR Lyrae stars or the non- 
variable stars bordering the instability strip — and the main-sequence turnoff, AV™, as 
illustrated in Figure 1. From the very first theoretical models of horizontal-branch stars 
having helium-burning cores and hydrogen-burning shells (Faulkner 1966), it has been 
known that the luminosity of the horizontal branch is only a very weak function of total 
stellar mass and hence of age, with the main effect of increasing total mass being to shift 
the star to the red along the zero-age horizontal branch (ZAHB). The mass of the helium 
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core, which has some influence on the luminosity of the ZAHB, is believed to be only a 
very weak function of the total stellar mass, and hence age. Therefore, the horizontal 
branch should be a good magnitude fiducial from which to measure the luminosity of the 
main-sequence turnoff, which fades at a rate of some 0.01 mag for each 1% increase in the 
age (e.g., Green, Demarque, and King 1987; Bergbusch and VandenBerg 1992). 

Surveys of GC ages using this vertical technique in a strictly differential sense - 
the so-called "AV" method — were carried out by Gratton (1985) and Peterson (1987), 
who came to conflicting conclusions: Gratton found from a sample of 26 clusters that the 
age spread in the inner halo was negligible, ^1.5 Gyr, with a mean age of some 15 Gyr 
(i.e., an age range ^10%), while that in the outer halo was larger, at least 20-30%, in 
essential agreement with the SZ model. Peterson, on the other hand, concluded that the 
age spread in his sample of 41 clusters could not be distinguished from the observational 
noise inherent to the technique, which amounted to ±17% or so. Peterson attributed 
the differences between his conclusions and Gratton's to three causes. First, his own 
sample was based on a complete survey of the literature, while Gratton's was selected 
from a relatively small set of papers; as a result, different numerical values of AV^^f 
were adopted for some clusters in the two surveys. Second, Peterson adopted a slightly 
different metallicity scale for the most metal-rich clusters, noting, however, that by itself 
this difference would account for no more than 3% in age. Finally, Gratton had calibrated 
the absolute magnitudes of the horizontal branch as a function of metallicity in two ways: 
(a) on the basis of theoretical stellar models for HB stars (Caloi, Castellani, and Tornambe 
1978), dM v /d[Fe/H] = +0.16, and (b) on the basis of the observed period-luminosity- 
amplitude relation for cluster RR Lyraes (Sandage, Katem, and Sandage 1981; Sandage 
1982), dMy/d [Fe/H] = +0.315. Gratton expressed a preference for the latter because it 
produced a greater consistency among the ages for clusters in the inner halo (in a sense, 
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assuming that which he was trying to determine) . Peterson rejected the latter metallicity- 
magnitude relation because it had little independent observational or theoretical support, 
but noted that this difference tended to affect the absolute age scale rather than the 
variation of age with Galactocentric distance. In addition to the points raised by Peterson, 
it should also be remarked that Gratton did not consider the observational errors in the 
individual age determinations, most particularly that the data for the more remote outer- 
halo clusters might simply contain larger random errors than for the closer inner-halo 
clusters. Conversely, it is also possible that Peterson's more heterogeneous sample included 
some poorer data, thus helping the observational noise to swamp the signal in his study. 
However, Peterson did reexamine his conclusions on the basis of the 17 clusters with the 
best observations, and found no significant difference from the results based on his full 
sample. 

More recently, Sarajedini and King (1989) used the difference between HB and turnoff 
magnitudes to estimate the ages of 31 Galactic GCs. They considered three different 
metallicity-luminosity relations for the horizontal branch: (a) the "classical" constant- 
luminosity assumption, dMyjd [Fe/H] = 0; (b) a relation based on theoretical HB mod- 
els (Lee, Demarque, and Zinn, private communication — presumably those used by Lee 
et al. 1994 — with a zero-point adjustment to correct from Y = 0.23 to Y = 0.24), 
which was also roughly consistent with observed luminosities based on the Baade-Wesselink 
method, dMyfd [Fe/H] = +0.17; and (c) Sandage's (1982, and references therein) period- 
luminosity-amplitude relation, dMy/d [Fe/H] = +0.35. As with Gratton's (1985) study, 
both the absolute ages and the slope of the age-metallicity relation were found to be sensi- 
tive to the choice of metallicity-luminosity law for the horizontal branch, with luminosity 
independent of metal abundance implying a strong age-metallicity relation and Sandage's 
relation implying none. All HB relations implied an age distribution with a full-width at 
half-maximum of roughly 45%, implying a one-sigma dispersion of 19%, and a total range 
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in excess of 60%, resembling Zinn's (1993) ongoing accretion model even more closely than 
SZ's slow contraction model 2 . 

Finally, the AV method has very recently been applied to 43 Galactic globular clus- 
ters by Chaboyer, Demarque, and Sarajedini (1996; CDS). Adopting a range of slopes 
dMy/d [Fe/H], they derive age ranges for the Galactic globular cluster population. For a 
reasonable intermediate metallicity slope, dMy/d [Fe/H] = +0.20, they infer an age range 
(±2<t) of some 9 Gyr for the Galactic globular clusters as a whole, with a mean cluster age 
of some 18 Gyr. However, they recognize that four of the clusters in their sample (Ter 7, 
Pal 12, IC4499, and NGC6652) are considerably younger than the rest and may repre- 
sent a distinct population; removing these clusters from the sample, they derive a reduced 
age range (again ±2cr) of ~5 Gyr {i.e., a = 1.25 Gyr) for the remainder. To arrive at 
these numbers, they have subtracted (in quadrature) the estimated observational errors in 
[Fe/H] and AV^q from the observed dispersion in the apparent ages; these are estimated 

2 However, it is possible that the inferred full- width at half- maximum has been inflated 
somewhat by the authors' use of the "generalized histogram" to characterize the distribu- 
tion of cluster ages (their Fig. 1). The age distribution of the Galactic globular clusters 
has some true, intrinsic width. The observed distribution will be somewhat broader than 
the true one, due to random observational errors. In creating the generalized histogram 
the observed distribution is convolved with a Gaussian kernel having a standard deviation 
of the same order of magnitude as the observational errors. The generalized histogram 
has thus been broadened by the error distribution twice, once by Nature and once by the 
astronomer. A perhaps more valid estimate of the true intrinsic age dispersion would be 
obtained by subtracting the variance of the observational errors from the sample variance; 
in this case, the intrinsic age dispersion of the halo clusters would be reduced from 19% to 
8% (1.1 Gyr for a mean age of 14.3 Gyr) if Sarajedini and King have correctly estimated 
the observational errors. 
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to amount to ±1.6 Gyr. The observational uncertainties in AV^^f , in particular, were 
calibrated by comparing independent determinations of /SV^q for 13 clusters, from which 
CDS concluded that the actual external errors of the published magnitude differences were 
in fact only 0.61 times as large as the original authors had estimated. The conclusions 
of the paper rest heavily on this assertion that numerous independent investigations have 
systematically overestimated their observational errors by some 64%: the observed disper- 
sion in the cluster ages including the four anomalously young ones is only ~3 Gyr; it is the 
quadrature difference between this and the typical estimated age uncertainty of 1.6 Gyr 
that leads to an inferred value of ~2.25 Gyr for the intrinsic age dispersion in the halo. 
Had the uncertainties in AV^q given in the original papers been adopted instead, the 
uncertainty in the individual ages would have been taken to be typically some 2.6 Gyr, 
implying an age dispersion of a/ (3) 2 — (2.6) 2 ~ 1.5 Gyr including the four anomalously 
young clusters, and allowing for no intrinsic age dispersion at all if those four clusters are 
set aside. 

And it is not clear whether this "correction factor" — the ratio of perceived external 
to internal error — is well justified. It is based on multiple studies of a few well-observed, 
fiducial clusters that the authors have tabulated and discussed in an Appendix. Starting 
at the top of their list, three citations are given for 47 Tucanae: Hesser, et al. (1987); 
Chaboyer, Sarajedini, and Demarque (1992); and Sarajedini and King (1989; SK). In 
going back to those sources, we find that the AV^q published by Chaboyer, Sarajedini, 
and Demarque is taken from Buonanno, Corsi, and Fusi Pecci (1989; BCF) who in turn 
state that their value is based on the photographic data of Lee (1977), while SK state that 
their value for AV^^ is itself based on the data of Hesser, et al. (1987). The second 
cluster on their list, NGC 288, has citations to BCF, SK, and Bergbusch (1993) taken in 
conjunction with Pound, Janes, and Heasley (1987). Again going back to these sources, 
we find that BCF cites previous photographic work by Buonanno and collaborators for 
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the horizontal-branch photometry, while the turnoff magnitude was determined from their 
own CCD data. SK list Olszewski, Canterna, and Harris (1984) and — again — Pound 
et al. (1987) as the source of the data used in their analysis. It is difficult to understand 
how the mixing of photographic and CCD data, and how the web of interlocking citations 
will have affected the validity of this "correction factor" of 0.61. Grundahl (1996) has 
compared the "vertical" ages of CDS with the "horizontal" (see below) ages of Richer et 
al. (1996). He concluded, "This comparison revealed large discrepancies between the two 
methods and the likely source found to be [sic] due to errors in the adopted AV^^f values 
from the literature ... It was also argued that this indicates that different authors do not 
determine the location of the TO in consistent ways and that one is most easily led to an 
underestimation of the true values." 

Whatever one concludes about the reality of the intrinsic dispersion of globular cluster 
ages, CDS found, as did Gratton and Sarajedini and King before them, that for small 
values of dMv/d [Fe/H], there is a distinct tendency for more metal-rich clusters to be 
younger than metal-poor ones. CDS quantify this as there being a high probability (> 
98% confidence level) of metal-rich clusters being younger than metal-poor ones for all 
values dMv/d [Fe/H] < 0.24 (their Fig. 2); however, according to them, this trend is 
superimposed upon a large apparent age spread at each metallicity. No compelling evidence 
was found for a trend of age with Galactocentric distance, although it must be noted that 
their sample did contain only one cluster outside a Galactocentric radius of 22 kpc. 

The "horizontal" method of measuring ages offers an alternative to the vertical tech- 
nique. VBS and Sarajedini and Demarque (1990; SD) formalized some folk wisdom that 
had been understood by the stellar evolution community for decades: within limits, stellar 
evolution tracks and model cluster isochrones are rather self-similar. Differences in mass 
and chemical abundance produce loci that differ in their absolute vertical and horizontal 
position in the theoretical Hertzsprung-Russell diagram but, in fact, their shapes do not 
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differ by very much. By shifting theoretical isochrones both vertically and horizontally in 
the H-R diagram in order to register them at the magnitude and color of the turnoff, VBS 
demonstrated graphically that the principal morphological difference between isochrones 
was a shrinking of the color difference between the turnoff and the lower part of the giant 
branch with increasing age (see their Fig. 2). At least for the lower- metallicity clusters 
([Fe/H] < —1.2 or so), the effect of metal abundance appeared to be negligible (see VBS's 
Fig. 3) , although the authors did stress that it would be dangerous to apply this technique 
to the comparison of clusters with greatly different metallicities, because the errors in the 
stellar models — due, e.g., to our imperfect understanding of convection and opacities 
- might themselves depend on Z. Therefore, given a sample of clusters of similar metal 
abundance, by simply measuring the differences in color between the bluest point at the 
turnoff and the lower part of the giant branch in some self-consistent way, one has a direct 
measure of the spread in age within the sample. 

VBS demonstrated that the most metal-poor of the nearby clusters with the best pho- 
tometry, NGC4590 = M68, NGC6341 = M92, NGC6397, and NGC 7099 = M30, with 
[Fe/H] ~ —2, show no perceptible range in age (Figure 2): the measured dispersion corre- 
sponds to ~0.3 Gyr, which is consistent with the measuring errors. Furthermore, consider- 
ing this group of GCs, stellar samples taken from different parts of the same cluster show 
variations like those found between different clusters, again suggesting that these differ- 
ences are consistent with the measuring errors. Subsequently, Heasley and Christian (1991) 
and Durrell and Harris (1993) have reported that the ages of NGC 5024, NGC 5053 and 
NGC 7078 = M 15, too, are indistinguishable from that of these four clusters. Conversely, 
among the most-metal rich of the clusters that they considered ([Fe/H] ~ -1-3), VBS 
concluded that NGC 288 was probably some 2.5±1.5 Gyr older than its second-parameter 
counterpart, NGC 362, and most of the other clusters in that abundance range (NGC 1261, 
NGC 2808, NGC 3201, NGC 5904 = M5, and Palomar 5). The famous anomalous cluster 



19 

Palomar 12 may lie just outside the abundance range where the horizontal method may 
be applied with confidence, but when the method is used nonetheless, an age difference 
of some 25% between Palomar 12 and more typical clusters is confirmed. At intermedi- 
ate metallicities ([Fe/H] ~ -1.6), the situation was rather less clear, due primarily to a 
shortage of good published data. Among most of the clusters in this sample (NGC 5272 
= M3, NGC 6218 = M12, NGC 6752, and NGC 7492), there was no compelling evidence 
for a range of age. The comparison of M3 with NGC 6752, in particular, suggested that 
age is not responsible for the notable difference in the horizontal-branch morphologies of 
these two clusters. Conversely, a comparison of M 3 with some of Stetson's unpublished 
data for NGC 6205 = M 13 suggested that the latter may be older than the former by some 
1-2 Gyr, which is of the correct sign but falling somewhat short of the 3 Gyr difference 
predicted from their HB morphologies. So the conclusions of VBS were: first, that the 
age spread in the available sample appeared to be unmeasureably small among the most 
metal-poor clusters: almost certainly ^0.5 Gyr, or some 3%, which would be a remark- 
able coincidence if the true duration of the cluster-formation epoch was of order 5 Gyr, 
or ~40%. Second, among the intermediate-metallicity clusters, there was a hint of the 
possibility — no more — that age might (sometimes?) be a second parameter; i.e., the 
results were ambiguous. Third, among the more metal-rich clusters, there was stronger 
evidence for a range of age, up to as much as 25-30% if Pal 12 represents one extreme 
of the age range for normal clusters, or as little as 2-3 Gyr (12-20%) if Pal 12 is a true 
anomaly with an origin quite different from that of the main halo cluster population. SD 
applied their own version of the horizontal method specifically to the NGC 288-NGC 362 
comparison, and likewise concluded that NGC 288 is older (by some 2.5-3 Gyr). 

Each of these approaches, the vertical and the horizontal, is doubly differential: the 
difference between two things measured in one cluster is compared to the same difference 
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measured in another. Both are independent of distance, reddening, and errors in the as- 
sumed photometric zero points. Apart from this common strength, the advantages and 
disadvantages of the two methods are different. The vertical method is on a sounder theo- 
retical footing, because the luminosity of the turnoff depends primarily on properties and 
processes in the stellar core, where the physics is simple and comparatively well under- 
stood. By contrast, the horizontal method is affected by the stellar radii, and hence might 
be sensitive to less-well understood envelope physics such as opacities and convection: if 
these were to differ from our understanding of them in ways which vary from cluster to 
cluster, then our age estimates would be affected. However, Figure 3 shows the same data 
as are plotted in Fig. 2, but compares them to a set of isochrones which are identical 
except for the choice of mixing- length parameter: ckmlt = 2.5, as opposed to 1.5. The 
comparison of Fig. 3 with Fig. 2 demonstrates yet again that, while the horizontal method 
cannot at present be used to determine absolute ages, a systematic error in our knowledge 
of convection theory at some particular metal abundance will not endanger conclusions 
about the relative age range at that abundance, unless there is cluster-to-cluster variation 
in some other parameter that alters the convection mechanism. Similarly, Figure 4 shows 
that the horizontal method is not particularly sensitive to unrecognized cluster-to-cluster 
differences in overall metallicity (for [Fe/H] < -1.2 or so), oxygen- and other a-element 
relative abundances, or helium abundance, unless a secondary effect of these hypothetical 
differences is to alter the physics of the envelope convection. But, as discussed in Vanden- 
Berg et al. (1996), it has become apparent that perceived turnoff magnitudes may also be 
subject to the vagaries of modelling stellar T eJ j's, and hence the vertical method may not 
be completely immune to these hidden-parameter uncertainties either. 

However, all such theoretical concerns aside, the apparent magnitude of the bluest 
point at the turnoff can be quite hard to measure, because the stellar locus is (by definition) 
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vertical at that point. Theoretical models predict that the color of the main sequence in B 
V will be constant to <0.01 mag over a range of ~ 0.6 mag, and constant to <0.02 mag over 
a range of nearly 0.9 mag. (VBS have shown that these numbers do not depend sensitively 
on age or metal abundance.) Thus, slight calibration errors, stochastic arrangements of 
stars in the diagram, the inclusion of composite-light binaries or marginal blue stragglers 
in the sample, or the preconceptions of the astronomer hand-drawing the mean locus, can 
easily bias the estimate of the vertical position of the bluest point by a meaningful amount. 
Conversely, the verticality of the turnoff makes it exceedingly easy to estimate the color 
of the bluest point. By effecting the vertical registration of two cluster sequences not by 
the apparent magnitude of that bluest point, but by the magnitude of the point on the 
main sequence at some fixed A (color) redder than the bluest point, where the slope of 
the cluster locus is appreciable, the VBS approach greatly reduces the uncertainty in the 
vertical offset which, since the turnoff is vertical and the lower RGB is nearly so, is of 
significantly reduced importance in any case. 

The vertical method relies on a fine measurement of the apparent magnitudes of 
objects that differ in luminosity by 3-3.5 mag (cf. Fig. 1), a factor of order 20. Early on, 
many observing programs were designed to do one of two things, but not both: either to 
define the horizontal branch and the giant branch down to some level below the HB, in 
which case the turnoff stars were undetected, or to measure the turnoff and as far down 
the main sequence as possible, in which case fields were chosen with a minimum number 
of giants and HB stars, and/or those brighter stars were saturated. Many of the earlier 
values of AV therefore were based on combining two different data sets, with the shallower 
data often coming from photographic work. Programs intended to bracket the HB and the 
turnoff within a single data set have more recently been undertaken, but it still requires 
a major effort to obtain tractable samples of both turnoff and horizontal-branch stars in 
the same images, and such studies are still vulnerable to non-linearities in the detector 
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and to possible subtle differences in the way the data analysis software treats bright, 
well-exposed stars as compared to faint, noisy ones. (For instance, in synthetic aperture 
photometry, bright stars are optimally measured in large apertures, while faint stars are 
better measured in small apertures, with an empirical correction to large-aperture values 
(e.g., Stetson 1989). In profile-fitting photometry, the dominant source of uncertainty in 
the measurement of bright stars is point-spread- function mismatch; the measurement error 
of faint stars is dominated by photon noise and the estimation of the local sky brightness.) 
In the horizontal method, the turnoff stars are compared to stars that are both much more 
like them in apparent magnitude (cf. . Fig. 1) and considerably more common than HB 
stars. Furthermore, the comparison is in color and any residual systematic nonlinearities 
in either detector or analysis will tend to cancel when the magnitudes are differenced to 
form the instrumental color index. 

The purely observational uncertainties in the horizontal method can be reduced sim- 
ply by using colors with a longer wavelength baseline (Stetson 1993). For instance, B I is 
roughly 2.5 times more sensitive to temperature than B-V, which produces the straight- 
forward factor of 2.5 improvement on the leverage for estimating temperature differences 
given fixed photometric errors by, in effect, increasing the curvature of the turnoff by a 
factor 2.5, and increasing the relative displacement of the lower RGB from the turnoff by 
the same factor. The use of B I also greatly aids in distinguishing the true main sequence 
from binaries and marginal blue stragglers. Furthermore, in B I the main sequence is more 
horizontal and the giant branch is more vertical than in B-V, thus rendering the exact 
vertical registration of two cluster sequences even easier and even less important. There is 
no analogous way to increase the power of the vertical method, since the turnoff and the 
flat part of the horizontal branch have similar temperatures, so the magnitude difference 
is roughly constant at all wavelengths. The similar colors of turnoff and HB stars free 
the vertical method from errors in the color transformation from instrumetal to standard 
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photometric indices. Still, with a color difference of order 0.25 mag (in B-V) between the 
turnoff and the base of the giant branch, the horizontal method is not overly sensitive to 
transformation errors of a percent or so, particularly when used to intercompare clusters 
observed during the same observing run. 

A final point about the vertical method is that it breaks down for ages that are too 
small. This is because, for horizontal-branch masses greater than 0.9-1.0 Mq, the ZAHB 
turns upward and then returns toward the blue, brightening all the while {e.g., Faulkner 
1966, Demarque and Hirshfeld 1975). The correctness of the models in this prediction is 
dramatically illustrated by the CMD of the dwarf spheroidal galaxy in Carina published 
by Smecker-Hane et al. (1994; their Figs. 2 and 3), where the clump of core helium- 
burning stars belonging to a ~6 Gyr population lies roughly a quarter-magnitude above 
the horizontal branch of a ~15 Gyr population. Therefore, below some age (which is greater 
than 6 Gyr), the vertical method becomes ambiguous, as a given AV can correspond to 
two possible ages. There is no corresponding barrier to the use of the horizontal method 
at younger ages. On the other hand, VBS stress that their horizontal method cannot 
be applied with confidence to clusters more metal-rich than some minimum value lying 
in the neighborhood of —1.2: for more metal-rich clusters the color difference between 
the turnoff and the base of the giant branch depends upon both metallicity and age. 
Perhaps the horizontal method can still be applied to determine differential ages for metal- 
rich clusters shown spectroscopically to have closely similar abundances. Alternatively, 
Buonanno et al. (1993) have demonstrated that enforcing consistency between vertical 
and horizontal techniques can simultaneously constrain age and metallicity differences, at 
least toward the metal-poor end of the globular-cluster abundance distribution; perhaps 
the same sort of approach can be made to work for more metal-rich clusters. However, 
the propagation of accumulated random errors may limit the ultimate precision achievable 
with such compound, multiply differential techniques. 
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Other types of vertical approaches are possible. For instance, given the theoretical 
prediction that the luminosity at the helium flash is nearly independent of the mass, and 
hence the age, of a star — certainly for any mass within the range 0.7 < M/ Mq < 1.0 (cf. 
Sweigart and Gross 1978) — the magnitude difference between this reference luminosity 
and that of the turnoff can be used to determine relative cluster ages. Of course, the 
brightest giant in a globular cluster will not be precisely at the RGB tip, but Monte Carlo 
simulations by, e.g., Crocker and Rood (1984) have shown that with as few as 45 stars in the 
upper 2.5 mag of the giant branch, the brightest should be within 0.11 mag of the helium 
flash luminosity 68% of the time. (For this reason, the helium-flash luminosity cannot be 
used as a reference in the Galactic globular clusters of the lowest mass, because some of 
them do not contain as many as 45 giants.) Making use of the extensive surveys of the 
bright giants in a large sample of GCs that had already been carried out (Frogel, Persson, 
and Cohen 1981, 1983), VandenBerg and Durrell (1990) vertically shifted the CMDs for 
several globulars having similar metallicities until the brightest giant in each cluster had 
the same V magnitude, and then compared their respective turnoff magnitudes. Those for 
NGC288, NGC362, and M5 were found to be sufficiently similar that VandenBerg and 
Durrell concluded that the ages of these three cluster were "not detectably different" (i.e., 
any variation was less than the estimated 1-2 Gyr la uncertainty). The same conclusion 
was reached for the group of very metal-poor GCs comprising M 15, M 30, M 68, and M 92. 
(It is encouraging that these results are close to those obtained subsequently by VBS, who 
used a technique that is capable of much higher precision.) 

But this approach has many practical difficulties. For one, as already implied, the 
brightest giants are very rare, and the method is therefore sensitive to small number 
statistics. For another, the V magnitude of the RGB tip is a fairly strong function of the 
metal abundance (see, e.g., Fig. 1; also VandenBerg 1992), so that even relatively small 
differences in chemical composition could, if not taken into account, lead to misleading 
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results. It would be far better to use infrared photometry: for instance, Da Costa and 
Armandroff (1990) have proposed that the RGB-tip /-band luminosity is constant (or 
very nearly so) in GCs. In addition, giant-branch tip stars are far more luminous than 
even horizontal-branch stars and, furthermore, they are extremely red: the use of a long- 
wavelength filter (such as I) would further amplify the magnitude difference between the 
tip stars and the much bluer turnoff, greatly compounding the dynamic-range problem 
with the vertical method. 

Fusi Pecci et al. (1990), and subsequently Sarajedini and Lederman (1991), have 
proposed that the magnitude of the "bump" in the giant-branch luminosity function can 
be used to estimate GC distances. It would follow that the magnitude difference between 
the bump and the main-sequence turnoff could also be used as a distance-independent age 
indicator. The practical difficulties with this method, too, would be large. Apart from the 
fact that the observed magnitude of the bump is not yet successfully predicted by theory, 
in most cases it would be necessary to survey a major fraction of a given cluster to obtain a 
sample large enough to define the bump accurately. Much of this area would be too crowded 
to measure the turnoff stars with adequate precision, so the sample of "bump" stars and of 
turnoff stars would be measured under significantly different conditions. Even then, unlike 
the horizontal branch, the bump differs from the turnoff in color as well as magnitude, so 
the age index would be sensitive to erors in the color transformation. Anyway, with the 
current generation of detectors, these magnitude fiducials will more probably be used in 
the same way as Sandage (1953) used the horizontal branch, as a means to estimate the 
absolute distance to a cluster, from which the apparent magnitude of the main-sequence 
turnoff as determined by a different set of observations can be used to infer the age, an 
approach which includes the difficulty of making sure the two studies are truly calibrated 
to the same system and, hence, is not purely differential. 
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5. A TEST OF THE HYPOTHESIS THAT AGE IS THE SECOND PARAMETER 
In their latest paper, Lee, Demarque, and Zinn (1994) presented a strong case that age 
must be the cause of the observed variations in HB morphology among clusters of the same 
[Fe/H] (the "first parameter") and as a function of Galactocentric distance (Rc)- They did 
so by the process of elimination; that is, they demonstrated that variations in any of the 
other usual second-parameter suspects — helium abundance, CNO abundance, and core 
mass (Rood 1973, Renzini 1977) — lead to inconsistencies with observations. Furthermore, 
they argued that the recently discovered correlations of the length of the blue HB tail with 
cluster density and concentration (Fusi Pecci et al. 1992, Buonanno 1993), which might 
be suggesting that environmental influences drive differential mass loss, is not apparent 
in the observed trend of HB type with Rq- Only the age explanation appeared to be 
comparatively trouble-free, and in view of the indisputable fact that young GCs do exist 
in the Galaxy, the evidence seemed compelling to them that age is the dominant (possibly 
the sole) second parameter. 

Figure 5 shows a composite version of one of the diagrams contained in the Lee et al. 
(1994) paper (their Fig. 7). It plots the HB type, as quantified by (B -R)/(B + V + R) - 
where B, V, and R represent the numbers of blue HB stars, RR Lyrae variables, and red 
HB stars, respectively — as a function of [Fe/H] for 83 GCs. Different symbols are used to 
indicate the division of the clusters into three radial zones: the solid lines illustrate the Lee 
et al. synthetic HB calculations assuming that all stars undergo the same amount of mass 
loss during the giant-branch evolution and that only [Fe/H] varies from cluster to cluster. 
The three curves are effectively "HB isochrones," with the upper/lower ones showing how 
the position of the middle locus would be changed by an age increase/decrease of 2 Gyr. 
This is admittedly a very tantalizing figure, showing as it does how the spread in HB type 
at a given [Fe/H] and the variation in HB morphology with Rq can be accounted for by 
differences in age. 
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However, as already noted, there are some difficulties for the pure age hypothesis. 
Using their A(color) technique, VBS derived differences in age between NGC288 and 
NGC362 and between M3 and M 13 that are smaller than the ones implied by Fig. 5. 
According to the "HB isochrones" computed by Catelan and de Freitas Pacheco (1993), 
for any reasonable estimate of the absolute cluster ages, NGC 288 must be ^ 5-6 Gyr older 
than NGC 362, if age alone is responsible for the differences in their HB morphologies. Lee 
(1991) has claimed that he could reduce the age difference to ~3-4 Gyr if the stars in both 
clusters had [O/Fe] = 0.85, but such a high oxygen abundance has pretty well been ruled 
out observationally: Dickens et al. (1991) have obtained [O/Fe]^+0.2 in unmixed stars in 
both NGC 288 and NGC 362; Sneden et al. (1991, 1992, 1994) and Kraft et al. (1992) have 
obtained similar values in clusters both more metal-rich and more metal-poor than these; 
Bell, Briley, and Norris (1992) have argued that CO band strengths in the metal-poor GC 
NGC 6397 are incompatible with O abundances as high as [O/Fe] = +0.6; and Brown and 
Wallerstein suggest that [O/Fe] « +0.5 over the entire range of iron abundances found 
in lo Centauri. Lee has also remarked that a large difference in [Fe/H] between the two 
clusters could accomplish the same thing, but this suggestion appears to be untenable as 
well (see, e.g., the spectroscopic studies by Caldwell and Dickens 1988 and Dickens et al. 
1991). In addition, based on the AV method, Stetson et al. (1989) have inferred that 
NGC 362 and M 5 have very similar ages: the same conclusion, using the same approach, 
was reached by Catelan and de Freitas Pacheco (1995) for the M 3-M 13 pair. These results 
clearly conflict with the expectations from Fig. 5, as does the recent work by Stetson and 
VandenBerg (1996) on M2 and M3. The latter have found, from differential photometry 
obtained on the same night with the same telescope, that the main-sequence to lower- 
RGB CMDs of these two clusters, which also form a second-parameter pair (see Fig. 5), 
are virtually identical, implying indistinguishable ages. 
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There are at least three further problems for the age interpretation of the variations in 
HB type among GCs. First, taken at face value, Fig. 5 suggests that all of the clusters with 
extremely blue HBs, including especially the intermediate-metallicity cluster NGC 288, are 
much older than all of the extremely metal-deficient systems like M 15 and M92. This 
seems highly improbable, as also noted by Bergbusch (1993), who found that isochrone 
fits favored very similar ages for NGC 288 and M 92. Second, among the "young" globular 
clusters, Arp 2 has a "normal-looking" blue HB for its metallicity, in spite of its relative 
youth (see Buonanno et al. 1994). The same comment could be made concerning IC 4499 
(Ferraro et al. 1995), whose HB looks very similar to that in M3 and which has a similar 
metallicity, but is certainly younger. Some additional parameter must be at work which is 
compensating for the effects of younger age on the HB morphologies of these two clusters. 
And third, the pure age explanation cannot be reconciled with the observations of bimodal- 
HB clusters like NGC 1851 (Walker 1992) and NGC 2808 {e.g., Rood et al. 1993), which 
possess both very blue and very red HB populations with few stars between them. 

This brings us to our simple test of the age hypothesis: let us use the CMD for 
a bimodal-HB globular cluster as a template with which to compare those of NGC 288 
and NGC 362. To be specific, let us force a coincidence of the turnoff luminosities of all 
three clusters, to be consistent with the working hypothesis that they are coeval, and 
then intercompare their bright-star populations to check the validity of that assumption. 
We have chosen to use the NGC 1851 CMD as obtained by Walker (1992) to act as this 
bridge because its photometric sequences are much better defined than existing data for 
NGC 2808, which has a similar metallicity. 

Figure 6 shows the superposition of the upper main sequences, the subgiant branches, 
and the lower RGBs of NGC 288, NGC 362, and NGC 1851 when the horizontal and vertical 
offsets listed in the caption are adopted. There is clearly some uncertainty in the choice 
of these offsets, given the obvious differences in the cluster fiducials, but by making the 
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age-dependent subgiant branches coincide, we remove most (if not all) of the ambiguity 
inherent in the matching of turnoff luminosities alone. Figure 7 illustrates the upper 
part of the NGC 1851 CMD on which hand-drawn envelopes to the stellar distributions 
are indicated by the solid curves. Then, in Figure 8, these same hand-drawn boundaries 
are superimposed onto the bright-star populations of NGC 288 and NGC 362 using those 
offsets that were adopted in producing Fig. 6. 

The agreement is remarkable: NGC 1851 contains a red HB just like the one in 
NGC 362 and a blue HB very much like that in NGC 288. (Although the stars at the top 
end of the blue tail in NGC 288 appear to be somewhat brighter than those in NGC 1851, 
this could well be expected if the former evolved from fainter ZAHB locations than the 
latter, which would be consistent with the observation that NGC 288 has a longer blue 
tail than NGC 1851; cf. Fig. 1.) The main implication of this comparison is that all 
three clusters do, indeed, have the same age to within quite a small uncertainty (^ ±1 
Gyr). Moreover, small cluster-to-cluster differences in [Fe/H] or [ct/Fe] will not alter this 
conclusion because we have effectively used the Ay method, which is insensitive to mod- 
est changes in heavy-element abundances (see, e.g., Bencivenni et al. 1991, Caputo and 
Degl'Innocenti 1995). In addition, since NGC 1851 contains the NGC 288 and NGC 362 
CMDs, differences in Y or [CNO/Fe] cannot be the explanation of the diversity in HB 
morphology among GCs either, unless such differences occur among the stars of this one 
cluster (which seems unlikely). 

The NGC 288-NGC 362-NGC 1851 comparison is perhaps the strongest evidence so 
far - - and there have been many other indications (some noted above) - - that HB- 
morphology arguments alone (cf. SZ, Lee 1992, Lee et al. 1994) cannot be used with 
confidence to support a particular chronology for the formation of the Galaxy. It also 
shows how risky it is to apply the AV method to those systems that have no zero-age HB 
stars at the color of the instability strip. For instance, Sarajedini and King (1989) (also 
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see Chaboyer and Demarque 1994) inferred that NGC 288 was much older than NGC 362 
because, in their estimation, the magnitude difference between the HB and the turnoff was 
0.33 mag larger in the former cluster than in the latter. On the other hand, it is not cer- 
tain that the extension of the ZAHB appropriate for the red HB stars in NGC 1851 would 
coincide exactly with the ZAHB from which the blue HB stars evolved. For instance, if 
core rotation were more important in the bluer stars, then they may have started their 
core helium-burning phases with slightly higher core masses than those to the red of the 
instability strip. But the difference would be small: the reddest of the blue HB stars in 
NGC 288 must have evolved from much bluer and fainter ZAHB positions, as is suggested 
by the HB evolutionary tracks shown in Fig. 1 above. (A similar conclusion is drawn if 
comparison is made with NGC 2808, although these data are poorer.) 

But why have color-difference methods (VBS, SD) also found NGC 288 to be older 
than NGC 362? There may be a number of reasons. It is still possible that the two 
clusters do differ in age by ~l-2 Gyr given the inaccuracies in the available photometry 
(see Sandquist et al., 1996 for a recent comparison of NGC 288, NGC 362 and M5). Only 
with excellent photometry would the color-difference method be able to measure an age 
difference to this level of precision. In addition, the previous applications of the color- 
difference method to these GCs were limited by the small number of evolved stars that 
had been observed (a fact sometimes masked by the use of fiducial lines in comparison 
plots). Small-number statistics in conjunction with the known observational imprecision 
may well have dominated the A(color) measurements (recall that the formal error in 
the VBS determination of the age difference between NGC 288 andNGC362 was ±1.5 
Gyr). In this day of large- format CCDs, it would be relatively straightforward to improve 
the samples in both clusters by a large factor. Finally, NGC 288 is known to have a 
significant population of main-sequence binaries (Bolte 1992), as indicated by the strong 
redward asymmetry in the distribution of stars on the main sequence. This could cause the 
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estimated turnoff color to be redder than that of the single-star turnoff, thereby simulating 
an older cluster. 

We can only speculate about what else could be causing the variation of HB mor- 
phologies within the bimodal HB clusters, and from cluster-to-cluster at a given [Fe/H] and 
as a function of Rg- Core rotation may be an important factor, though it would be hard to 
explain the dependence of the second-parameter effect on position in the Galaxy in terms 
of angular momentum differences (see Norris 1981, 1983). The concentration of a cluster 
or its central stellar density may also play a key role in this phenomenon (Fusi Pecci et al. 
1992, and references therein). As noted by Renzini (1983), nearly all clusters with gaps 
on the HB have high central concentrations: close encounters may lead to the stripping of 
stellar envelopes or to the spinning-up of significant numbers of stars. Possibly increased 
angular momentum is a natural consequence of high concentration (Buonanno, Corsi, and 
Fusi Pecci 1985; Buonanno et al. 1986): according to Djorgovski and Meylan (1994), GCs 
closer to the Galactic center tend to be more concentrated with smaller, denser cores. But 
these ideas can hardly apply to NGC 288, which is a low-concentration GC, though gaps 
do seem to exist on its horizontal branch as well (see Fig. 8). However, perhaps the large 
binary population in this cluster (Bolte 1992, Bolte and Dubath 1996) is the critical clue 
to what is going on. The way in which the protoGalaxy collapsed may also have given rise 
to second-parameter effects (van den Bergh 1993). Indeed, it may well be that all of these 
contribute to a greater or lesser extent. 

The bimodal-HB clusters like NGC 1851 and NGC 2808 are unquestionably very im- 
portant Rosetta stones and they should be exhaustively studied. One certainly has the 
impression from Figs. 7 and 8 that whatever is afflicting NGC 1851 is also the cause of the 
difference in the HB morphologies of NGC 288 and NGC 362. 
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6. RECENT DEVELOPMENTS 
The launch and subsequent repair of the Hubble Space Telescope have finally permit- 
ted the determination of precise ages for globular clusters deep in the Galactic bulge and 
out at the extremes of the halo. Using the Planetary Camera of the original WF/PC on 
the unrepaired HST, Fullton et al. (1995) have obtained a color-magnitude diagram in the 
(V, I) bandpasses for the globular cluster NGC 6352 which, according to the dichotomy dis- 
covered by Mayall (1946) and Kinman (1959), and subsequently confirmed by Zinn (1985), 
belongs to the disk subsystem of globulars. By applying the AV method to NGC 6352 with 
respect to 47 Tucanae (also a cluster of the disk subsystem, but slightly more metal-poor 
and therefore requiring a small adjustment to AV to allow for the abundance difference), 
they concluded that NGC 6352 was the older of the two by 0.7±2.2 Gyr. Ortolani et al. 
(1995) present analysis of data obtained with WFPC2 of the near-Solar metallicity bulge 
clusters NGC 6528 and NGC 6553. Based on either the horizontal or vertical differential 
age indicators, these two clusters appear to be identical in age to ^ 10%. Comparing 
the clusters' luminosity functions to that measured for the bulge field in one of Baade's 
Windows, after registering the position of the HB, they find a very close correspondence 
of the cluster and field main sequences indicating nearly identical ages for these clusters 
and the field stars of the bulge. The authors also make a case that the clusters' ages are 
close to that of NGC 104 (47 Tucanae), which is down in metallicity by a factor of 4-5 
compared to NGC 6528 and NGC 6553. Because the effects of metallicity are very large at 
the metal-rich end of the scale, this conclusion is less easily justified. With that caveat in 
mind, it nevertheless appears that there is no evidence for a significant age range within 
this sample of objects from the inner Galaxy. Finally, recent work by Grundahl (1996) led 
to the conclusion that the moderate-metallicity thick-disk GC's M71, 47 Tuc, NGC 6352, 
and NGC 6760 have ages that are indistinguishable at the level of 1.5 Gyr, suggesting that 
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the oldest part of the thick disk formed over a very short time span. Conversely, the near- 
Solar metallicity clusters NGC 5927, NGC 6528, NGC 6553 formed some 3-4 Gyr later; 
following Burkert, Truran, and Hensler (1992), Grundahl suggests that these clusters may 
have formed along with the first generation of thin-disk stars. 

Stetson et al. (in preparation) have used the Wide Field and Planetary Camera 2 
on the post-repair HST to derive color-magnitude diagrams for the remote-halo clusters 
NGC 2419, Pal 4, and Pal 3. All three have the large core radii typical of the diffuse 
globular clusters of the outer halo, but while Pal 4 and Pal 3 both also have the very 
low masses characteristic of outer halo clusters, NGC 2419 is the third most luminous 
Galactic globular cluster known. In addition, while Pal 3 and Pal 4 share the anomalously 
red horizontal branch implying the operation of the second parameter, the horizontal 
branch of NGC 2419 is blue (normal for an inner-halo cluster with its measured metal 
abundance, [Fe/H] ps —2). The vertical and horizontal methods as applied to these clusters 
demonstrate that the age of NGC 2419 is indistinguishable at the level of 0.5-1 Gyr from 
that of NGC 6341 = M92, a cluster with similar metal abundance and HB morphology, 
while the age of Pal 3 is indistinguishable at a similar level of precision from that of M 5, a 
cluster with comparable metal abundance but having a very much bluer horizontal branch. 
The vertical and horizontal methods as applied to Pal 4 compared with M 5 are somewhat 
ambiguous, in the sense that both appear to suggest that Pal 4 is younger by several Gyr, 
but there is not a good correspondence of mean cluster loci at all magnitudes: whether this 
anomaly is real or due to some error in the photometric analysis is still under investigation. 
But the first hints from these studies are that (1) age is not the sole cause of the diversity 
in HB morphologies, and (2) the variation of mean age with Rq is small and possibly even 
zero. 



34 

7. SUMMARY 

As the result of the superb color-magnitude data that have been obtained through 
the use of CCD detectors on large telescopes and the Hubble Space Telescope, relative GC 
ages are starting to be much better known. The application of color- and magnitude- 
difference techniques to the observed CMDs is beginning to suggest that there is little or 
no dispersion in age, either at a given [Fe/H] or as a function of Galactocentric distance 
(Rg)- (There are undeniably a number of young GCs in the Galaxy, but they appear to 
represent a small and quite distinct minority and may well have been tidally pulled out of 
the Magellanic Clouds and/or a nearby dwarf spheroidal.) In this study, we have provided 
fairly strong (supporting) evidence that age is not the sole cause of the second-parameter 
phenomenon, if indeed it ever plays an important role at all. Consequently, the observed 
changes in horizontal-branch type with Rq can probably not be explained in terms of age 
variations alone. 

Whether or not GC ages vary with [Fe/H] is almost certainly the hardest of the 
relative-age questions to answer. To do so will require, among other things, the resolution of 
present distance-scale uncertainties (mainly concerning the variation of horizontal-branch 
luminosities with metal abundance), and an improved understanding of the chemistry of 
stars (notably the variation of [a/Fe], and especially [O/Fe], with [Fe/H]). We suspect that 
the confluence of various streams of thought, including but not limited to parallaxes and 
proper motions from HIPPARCOS, MACHO-spinoff studies of variable stars in the Large 
Magellanic Cloud and the Galactic bulge, HST observations of globular clusters in the 
Local Group, and advances in stellar evolution theory, will provide us with a trustworthy 
resolution of this question relatively soon. 

These matters aside, the recent results discussed in the previous sections suggest 
that the Eggen, Lynden-Bell, and Sandage (1962) picture of the formation of the Galaxy 
(especially the updated version described by Sandage 1990) seems to be a good description 
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of the dominant formation mode of the Galactic GC system that we observe today, but 
in addition there exists ample evidence that accretion/merger processes such as those 
envisioned by Searle and Zinn (1978) and Zinn (1993) have also contributed to the Galactic 
halo. We hasten to add, however, that no clear consensus has yet emerged — this remark 
represents only our interpretation of the available evidence as of mid-April, 1996. 
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FIGURE CAPTIONS 

Fig. 1. Bergbusch & VandenBerg (1992) isochrones for [Fe/H] = —2.03 and —1.26 and 
the ages indicated are plotted as solid curves from the bottom center to the upper 
right part of the figure, and the main sequence (MS), subgiant branch (SGB), and red- 
giant branch (RGB) evolutionary stages are labeled along the more metal-rich loci. 
For each metal abundance, the youngest age has the brightest, bluest turnoff and the 
bluest giant branch. Fully consistent zero-age horizontal branches (ZAHB) as given 
by Dorman (1992) are plotted as dashed curves, with the longer curve representing 
the ZAHB for the [Fe/H] = -2.03 case while, for the sake of clarity, only the red 
end of Dorman's ZAHB for [Fe/H] = —1.26 is shown (the lower of the two dashed 
loci). Solid curves rising from the metal-poor ZAHB locus represent a selection of 
horizontal-branch evolutionary tracks for different post-giant-branch stellar masses. 
The approximate location of the RR Lyrae instability strip is indicated by the parallel 
dotted lines. The meaning of the quantity AV^f — which represents the magnitude 
difference between the bluest point at the turnoff (TO) and the horizontal branch 
(HB) at the same color — is illustrated for the case [Fe/H] = -2.03. This quantity is 
clearly a function of age and somewhat insensitive to [Fe/H] given that the luminosity 
of the horizontal branch and of the turnoff point (at a fixed age) are both predicted 
to decrease as the metallicity increases (though not at exactly the same rate). 

Fig. 2. Panel (a): VandenBerg et al. (1996) isochrones for the specified parameters and 
ages in 2 Gyr increments from 12 to 18 Gyr (in the order indicated), after they have 
been registered so that they have the same turnoff colors and identical magnitudes 
at the point on their respective lower main sequences that is 0.05 mag redder than 
the turnoff. Panel (b): Registration of the M68, M30, and NGC 6397 fiducials to 
that of M 92 and their superposition onto the lower RGB segments of the isochrones 
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plotted in panel (a) . The sources of the data and the horizontal and vertical offsets 
that were applied to the cluster photometry are given by VandenBerg et al. (1990). 

Fig. 3. As for Fig. 2; in this case, the isochrones assume a value of «mlt = 2.5, but are 
otherwise identical to those portrayed in the previous figure. 

Fig. 4. A demonstration of how the shape of a 16 Gyr isochrone for [Fe/H] = —2.14, 
[a/Fe] = 0.3, Y = 0.235, and «mlt = 1-89 would change, according to VandenBerg 
et al. (1996) models, if the latter four parameters were varied by modest amounts 
(as indicated). The dashed and dot-dashed curves were registered to the solid curves 
at the magnitudes and colors indicated by the small arrows. 

Fig. 5. HB morphology vs. [Fe/H] for globular clusters having Rq < 8 kpc (closed 
circles), 8 < Rq < 40 kpc (open circles), and Rq > 40 kpc (crosses). The data and 
the solid lines, which represent theoretical isochrones produced by a synthetic HB 
code, are due to Lee et al. (1994). The clusters that have been identified on this 
plot are discussed in the text. 

Fig. 6. Comparison of the turnoff regions of the observed fiducial sequences for NGC 288 
(Bolte 1992), NGC 362 (VandenBerg et al. 1990), and NGC 1851 (Walker 1992b), 
after vertical offsets of +0.13 and —0.60 mag and color shifts of —0.018 and —0.035 
mag, respectively, were added to the tabulated sequences of the latter two clusters. 
Due to the obvious erratic difference among the sequence shapes, the vertical registra- 
tion was effected at the middle of the nearly horizontal part of the subgiant branch. 
As this part of the cluster locus is highly sensitive to age, any age differential among 
the clusters should cause a disregistration of their HBs when these offsets are applied. 

Fig. 7. Walker's (1992b) photometry for the red-giant and horizontal-branch stars in 
NGC 1851. Hand-drawn envelopes to the stellar distributions are delineated by the 
solid curves. 
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8. Overlay of the solid curves from the previous figure onto the Harris (1982) pho- 
tometry of the bright stars in NGC 362 and that by Bergbusch (1993) for NGC288, 
assuming exactly the same offsets as those adopted in producing Fig. 6. Zero-point 
shifts of —0.010 mag in B — V and +0.0635 mag in V were applied to the Bergbusch 
observations, given that he has found that these corrections are required to place his 
observations onto the system of the Bolte (1992) data, which were used to effect the 
registration shown in Fig. 6. At the red end of NGC 288's blue horizontal branch, 
stars appear to hug the upper envelope of the NGC 1851 horizontal branch; this may 
be because these stars have evolved from bluer zero-age horizontal-branch positions 
than the ones in NGC 1851, along tracks like those shown in Fig. 1. 
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